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ABSTRACT 


The shapes of liquid drops resting on flat surfaces were determined by a Runge- 
Kutta solution of the Laplace capillary equation. A characteristic length equal to the 
square root of surface tension divided by the product of density and gravity was used to 
nondimensionalize the numerical results. In addition, asymptotic solutions for small 
and large drops were combined to give explicit expressions for the maximum drop height 
and radius. These correlations apply for the complete range of liquid volumes and con- 
tact angles. 


LIQUID DROPS: NUMERICAL AND ASYMPTOTIC SOLUTIONS OF THEIR SHAPES 
by Kenneth J. Baumeister and Thomas D. Hamill 
Lewis Research Center 

SUMMARY 

The equilibrium shapes of liquid drops supported by a flat horizontal surface were 
determined as a function of the contact angle by solving the Laplace capillary equation. 
In addition, asymptotic solutions of the Laplace capillary equation for the case of large 
and small drops were combined to give explicit expressions for the maximum height and 
radius of the drop as a function of the contact angle for any liquid volume. 


INTRODUCTION 

The boiling and condensing heat -transfer characteristics of wetting and nonwetting 
liquids depends quite often on the geometrical characteristics of liquid drops in contact 
with solid surfaces. The contact areas are of particular importance in determining heat- 
transfer and vaporization rates. 

For example, in the two-phase-mist -flow regime in the liquid-metal Rankine boiler, 
spiral inserts centrifuge entrained liquid droplets to the heated-tube wall. At the wall, 
the shape of the drop must be known if the vaporization characteristics are to be pre- 
dicted. Poppendiek (ref. 1) has developed an analytical model for this flow based on a 
Leidenfrost boiling model. However, a complete analysis of Leidenfrost boiling re- 
quires a knowledge of the heat -transfer areas associated with a given size drop (ref. 2). 

These facts provide the incentive for the present report, in which the equilibrium 
shapes of liquid drops supported by a flat horizontal surface are determined. 

A liquid drop supported by a horizontal surface assumes a shape that depends on 
both molecular forces and body forces (gravity, etc. ). For example, the photographs in 
figure 1 (taken from ref. 3) show the effect of gravity on a mercury drop resting on a 
surface. Clearly, the gravitational body force g will affect the energy transfer to the 
drop in a heating situation. 


Molecular forces lead to the macroscopic phenomenon of surface tension. Since 
molecules are subjected to an unsymmetrical force field, the surface molecules are, in 
effect, pulled towards the interior of the liquid more strongly than towards the vapor. 
Complications arise because of the interaction between liquid molecules and the molecules 
of the solid supporting surfaces. These forces lead to the macroscopic phenomenon of 
contact angle, which is designated by 6 in figure 2. The effect of surface attraction is 
pictured in this figure. However, from a macroscopic viewpoint, the shape of the drop 
is determined by gravity, surface tension, liquid density, and contact angle. 

The mathematical formulation of the criterion of equilibrium of the drop surface is 
the Laplace capillary equation (ref. 4) 


1 ) 1 _ AP 

Rj R 2 <7 


(1) 


where or is the surface tension between the liquid and its vapor, and R^ are the 
principle radii of curvature at any point on the drop surface as shown in figures 3 and 4, 
and AP is the pressure difference across the liquid interface at the point in question. 
Equation (1) can be derived from static equilibrium considerations (ref. 5) or by a varia- 
tional technique of minimizing the total free energy of the drop (ref. 6). Equation (1) is 
rewritten in differential form by expressing the two radii of curvature in terms of differ- 
ential notation. 

Laplace, Gauss, Poisson, Gay-Lussac, and others (as cited in ref. 4) have at- 
tempted to solve equation (1). As yet, no general solution appears possible. The most 
complete work on this subject to date was by Bashforth and Adams (ref. 4). By suitable 
geometric transformations, Bashforth and Adams transformed equation (1) into a non- 
linear second-order ordinary differential equation. This differential equation was then 
solved numerically by hand and the results tabulated in reference 4. 

In their numerical solution of equation (1), Bashforth and Adams nondimensionalized 
the governing equation by introducing the radius of curvature of the drop at its apex as a 
characteristic length. However, for a given drop volume, the radius of curvature at the 
apex is not a known quantity. Consequently, as pointed out by Wark (ref. 7), parametric 
plots must be constructed from Bashforth and Adams tabulated results if the drop height 
and radius are to be explicitly related to the volume of a drop. This procedure proves 
to be inconvenient. 

The characteristic length 


X = 



( 2 ) 


2 



which is introduced in this report is more convenient to use since it depends only on the 
properties of the fluid and the gravitational field; whereas, Bashforth and Adams charac- 
teristic length depends on the shape of the drop. Because the characteristic length used 
in this report is different from that of reference 4, the numerical results were more 
easily obtained by a Runge-Kutta solution of equation (1) on a digital computer rather than 
attempting to use the tabulated results of reference 4. Also, the digital result has much 
greater accuracy. 

Numerical solutions are presented for the maximum drop radius, drop height, and 
contact radius in terms of the drop volume for contact angles from 1° to 180°. In addi- 
tion, asymptotic solutions of the Laplace capillary equation for the case of large and 
small drops are combined to give explicit expressions for the maximum height and radius 
of the drop as a function of the contact angle 9 for any liquid volume V. 


SYMBOLS 

C constant 

o 

g acceleration of gravity, cm/sec 

2 

g dimensional conversion factor, (g) (cm)/ (dynes) (sec ) 

(pertains only to work done in English system of units) 

h height of drop, cm 

h average height of drop 

h* dimensionless height of drop, h/X 

h* h/X 

h,^ maximum asymptotic (large drop) height of drop, cm 
h* dimensionless h^, h^/X 

j sin <p 

2 2 

AP pressure drop across interface, dynes/cm (N/cm ) 

R q radius of curvature of drop surface at apex of drop, cm 

R* dimensionless R Q , R Q /X 

Rj, Rg radii of curvature of drop surface, cm 

Rj, R| dimensionless radii of curvature, Rj/X and Rg/X 
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r 


maximum radius of drop, cm 
r* dimensionless r, r/X 

r c surface contact radius, cm 

r* dimensionless r , r /X 

C L C 

r^ radius of drop associated with h^, cm 
r* dimensionless r^, r^/X 

5 unit of arc length, cm 

S* dimensionless S, S/X 

3 

V volume of drop, cm 

V* dimensionless V, V/X** 

x radial coordinate 

x* dimensionless x, x/X 

Ax numerical integration increment, cm 
Ax* dimensionless Ax, Ax/X 

z height coordinate, cm 

z* dimensionless z, z/X 

z Q z at apex of drop, cm 

z* dimensionless z . z/X 

o o’ o' 

6 contact angle, deg 

X characteristic length (see eq. (2)) 

3 

liquid density, g/cm 

3 

P v vapor density, g/cm 

a surface tension, dynes/cm (N/cm) 

( p angle, deg (see fig. 4) 

Subscripts: 

i numerical integration index 

max maximum 
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GOVERNING EQUATIONS 


The solution of the Laplace capillary equation (eq. (1)) will determine the shape of 
the drop. This equation represents a force balance which relates the surface-tension 
forces and the pressure difference across the surface of a drop. At equilibrium the 
surface-tension forces must be balanced by the pressure difference across the free sur- 
face of the drop. 

By expressing the pressure difference across the drop surface in terms of a liquid 
head and by introducing the characteristic length \ into equation (1), the governing equa- 
tion becomes dimensionless and is of the form 



( 3 ) 


(The starred quantities denote dimensionless terms. ) The derivation of equation (3) is 
given in appendix A. 

Equation (3) must, of course, be rewritten with the radii of curvature expressed in 
differential form before it can be solved directly. If R| and R| are expressed in 
terms of x* and z*, equation (3) can be transformed into a second-order nonlinear dif- 
ferential equation (see appendix A) of the form 



( 4 ) 


where x* and z* are the dimensionless coordinates of a point P on the drop surface 
(see fig. 4). 

In principle, this differential equation must be solved to determine the shape of 
the drop. Since this equation cannot be solved explicitly, numerical procedures must be 
used. However, rather than a direct numerical attack on equation (4), reference 8 sug- 
gests the introduction of a new variable 


j = sin (p (5) 

as a means of simplifying the numerical work. When the new variable j is used, the 
Laplace capillary equation is transformed into a coupled pair of first-order differential 
equations (see appendix A). These equations are given by 
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( 6 ) 


dj = z* L 

dx* x* 


dz* 

dx* 



(7) 


Equations (6) and (7) represent the governing differential equations for the drop. 
These equations are solved by the Runge-Kutta method. (See appendix B for details of 
the numerical solution. ) 


DISCUSSION OF RESULTS 

The numerical Runge solutions of the Laplace capillary equation are presented 
graphically in figures 5 to 9. Basically, for a given dimensionless drop volume and con- 
tact angle, these graphs enable the calculation of the following parameters: 

(1) Maximum dimensionless height, h* 

(2) Maximum dimensionless drop radius, r* 

(3) Dimensionless contact radius, r* 

These three parameters, illustrated in figure 3, enable a sketch of the drop to be drawn 
and the contact areas to be determined. The numerical results from which these curves 
were plotted are given in table I. 


Drop Height 

In figures 5 and 6 the drop height h* and , respectively, are plotted as a 

function of drop volume for various contact angles from 1° to 180°. 

Any particular curve is characterized by two asymptotic regions. For very large 
volumes the drop height approaches a constant, the particular value of which depends on 
the contact angle. The second asymptote occurs for very small volumes. A very small 
drop assumes a spherical shape subject to the constraint imposed by a fixed contact angle. 
A completely nonwetting liquid (6 = 180°) forms a spherical drop in the limit as the 
volume approaches zero, while wetting liquids form truncated spheres in the zero 
volume limit. The two asymptotic regions are discussed quantitatively in the following 
paragraphs. 

In the large-volume regime, a unique relation exists between the members of the 
family of curves. That is, the asymptotic height is given in terms of the contact angle 
by (ref. 9) 


6 


h 


( 8 ) 


h* " 

y ag 

P L g 


= V 2(1 - cos 9 ) 


This equation was derived by neglecting one of the radii of curvature in the Laplace 
equation. Equation (8) also predicts that the slope in the large-volume regime should be 
-3, since 


or 


h* _ V2(l - cos 6) 
y*l/3 y*l/3 



( 10 ) 


The equation allows the asymptotic thickness of a drop to be computed for arbitrary con- 
tact angle in the large-volume regime. As seen in figure 5, the large-volume regime is 
determined by V* > 100. 

The vertical asymptotes of figure 6 characterize the small-volume regime. Each 
contact angle gives a different asymptotic value for the dimensionless drop height, which 
can be computed from elementary geometric considerations. The equation for the dimen- 
sionless height in terms of the contact angle, derived in appendix C, is given by 


lim 

V*-0 


h* 

v *l/3 


3(1 - cos 0) 


ni/3 


?r(2 + cos 9 ) 


(ID 


The small-volume asymptotes derived numerically for various contact angles in fig- 
ure 6 are identical to those predicted by equation (11), providing a further check on the 
accuracy of the numerical method. The plotted results of figure 6 indicate that a drop is 
"small" if V* < 0.01, which serves as the criterion for the applicability of equation (11). 

Drops whose volumes fall into the range 0. 01 < V* < 100 are neither spherical nor 
pancaked in shape. No theoretical expression for the height of drops exists in this inter- 
mediate range. Based on the two asymptotic solutions, however, an empirical expres- 
sion was deduced to fit the curves over the entire dimensionless volume range. The gen- 
eral correlation is given by 
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( 12 ) 


h* 


. 

4tt( 2 + cos 0)sin - 
2 

> . 
m 

V 

i 

co 

< 

* 

1 

] 


-l/3m 


where m is an empirical constant chosen to fit the data in the intermediate range. The 
best choice of m is m = 1. 0; this gave agreement to within 10 percent. Thus, 



3V* 


- 11/3 


3V* + 47t(2 + cos 0)sin - 

2 


(13) 


The general correlation is compared graphically with exact numerical results in figure 7. 
The particular form of equation (12) was chosen, because irrespective of the choice of m 
this equation approaches equations (8) and (11) in the limit of very large and very small 
volumes, respectively. 


Maximum Drop Radius 

In figure 8 the dimensionless drop volume V* is plotted against the dimensionless 
maximum radius with contact angle 6 as a parameter. These curves are also 

characterized by two asymptotic regions. 

For large volumes. 


V* = 7ir* 2 h* (14) 

leading to a slope of 2 on a log-log paper, since h* remains constant in this range. 
Solving for r* in terms of volume and contact angle by employing the expression for 
h* (eq. (8)) gives 


r 


* 

OO 


r 

* 

> 

l 

!/2 / a 

/ y* \ 

uV 2(1 - cos 0) 

v**j 


1/2 


( 15 ) 


For very small volumes, geometrical considerations (see appendix C) lead to 
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(16) 


= 77 jj(l - COS 0) 2 - - (1 - COS 0) 3 ^ 


R 


*3 


This equation gives a slope of 3 on log-log paper, as indicated by figure 8. 
Solving for the radius of the sphere gives 


R* = 


r *l/3 


* 1/3 [(1 - 


cos oy 


-(1 


cos 9) 


i 


1/3 


(17) 


For contact angles greater than 90°, the radius of the sphere represents the maximum 
radius of the drop. However, for contact angles less than 90°, the maximum drop radius 
r is given by 


r* = R* sin 6 = 


(cl/3 


sin 9 


C 1/3 [(l - 


cos 0) 2 - A (1 - cos 0) 3 


I 


1/3 


(18) 


Now, the large- and small-drop results can be combined in the intermediate region be- 
tween the two. 

In the intermediate range there is a gradual change of slope in which the two asymp- 
totic solutions blend into each other. There is no analytical expression for computing the 
maximum radius in the intermediate range. However, combining the two asymptotic so- 
lutions over the entire dimensionless volume range in the following manner gives, for 
9 < 90°, 




r 1//3 h* sin 2 0 


9 1 9 2/3 

cos By — (1 - cos 9y 
3 



l/2n 


and for 9 > 90°, 





,V3* 


(1 - cos 0) 2 


- - (1 - cos 9) 
3 


3" 



(19) 


( 20 ) 
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A choice of n = 1. 8 correlates the exact numerical results to within 2 percent. The 
agreement between the correlation and the exact numerical results is so close that no 
graphical comparison was made. 

As before, the particular choice of the correlating equation is such that independent 
of the value of n the equation will converge to the proper asymptotic solution for the 
case of large and small drops. 


Contact Radius 

The contact radius r*, which determines the area of physical contact between liquid 
and solid, is plotted against the maximum radius r* in figure 9 for contact angles of 
135° and 180°. For contact angles equal to or less than 90° the maximum radius equals 
the contact radius. 


Average Drop Thickness 

A parameter that has had some use in heat -transfer analysis is the average drop 
thickness h defined as 


h* 



and 


h* _ y* 2 /3 
y* l/ 2 j?r*2 

Both these quantities are listed in table I for various size drops. 


( 21 ) 


( 22 ) 


Experimental Data 

Some experimental data were available for various liquids, as indicated in figures 6 
and 8, on the relation between volume and the maximum drop radius and height. The 
data given in reference 10 are for liquids in Leidenfrost film boiling, in which the liquid 
does not wet the surface. Consequently, the contact angle is 180°. 
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CONCLUDING REMARKS 


Numerical solutions of the Laplace capillary equation were presented graphically for 
various contact angles. In addition, the generalized correlating equations for the maxi- 
mum thickness and radius of a drop for arbitrary drop volume and contact angle were 
synthesized by combining the asymptotic solutions for very large and very small drops. 
These correlation equations are accurate to within 10 percent on drop height and 2 per- 
cent on maximum drop radius. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, May 22, 1968, 

129-01-11-02-22. 
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APPENDIX A 


TRANSFORMATION OF LAPLACE'S CAPILLARY EQUATION 


The governing equation for the shape of the drop, Laplace's equation, is given by 
equation (1) (repeated here for convenience). 


_1_ + _1_ _ AP 
Rj R 2 cr 


(Al) 


The pressure difference across the drop is proportional to the sum of the reciprocal 
of the radii of curvature in two mutually perpendicular directions, as depicted in figure 4. 
Note in equation (Al) that the expression 



is an invariant with respect to the coordinate system chosen. 

The pressure difference AP in equation (Al) is a function of z which can be con- 
veniently expressed as 


AP = p L (z - z Q ) 1. + C 


(A2) 


where the vapor density p y is assumed much less than p^ and the constant C repre- 
sents the jump in pressure across the liquid-liquid interface at z = z Q and x = 0 (the top 
of the drop). 

At the top of the drop, the radii of curvature are equal; thus, Laplace's capillary 
equation (eq. (Al)) becomes 


AP 


z z o 
x=0 



C 


(A3) 


where R Q is the radius of curvature at the top of the drop. Substituting equation (A3) 
into the expression for AP(e.g., eq. (A2)) gives 


AP 



2 a 





(A4) 


12 



Up to this time, the quantity z Q represented an arbitrary distance from the refer- 
ence coordinate to the top of the drop. However, z Q is now taken as 


z o 


2 a 


R o 


(A5) 


Substituting this expression into equation (A4) gives 


AP = p T z -i- (A6) 

^ g 
5 c 

Thus, the governing equation (eq. (Al)) becomes, by substituting equation (A6) into equa- 
tion (Al), 



R 1 **2 


P L z 


<7 


(A7) 


The radius of curvature at the top of the drop R Q can now be written from equa- 
tion (A5) as 


R = ?£_ (A8) 

P L*of- 

®c 

The parameter z q is considered an independent variable which can be chosen arbitrarily 
to give various radii of curvature at the top of the drop. 

Choosing z Q as the independent variable, the volume of the drop becomes a depen- 
dent variable, that is 


V = f(z Q ) (A9) 

The functional form in equation (A9) cannot be determined explicitly, rather it can be 
represented in a geometrical fashion from the numerical solution of equation (Al). 

The governing equations have been discussed briefly and now will be transformed to 
a more convenient form for solution. 

Dividing both sides of the governing equation (A7) by the characteristic length X 
(defined by eq. (2)) yields 
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where 


(A10) 


(All) 


(A12) 

(The starred symbols represent dimensionless coordinates. ) 

The dimensionless Laplace capillary equation given by equation (A10) can now be 
written in terms of dimensionless variables x* and the angle <p (see fig. 4). 

Consider the point P on the surface of the drop defined by the coordinates <p and 
x*. The radii of curvature R| at the general point P can now be represented from 
geometrical considerations (see fig. 4) 

R* = _^!L = (A13) 

sin cp j 

The radius of curvature R| must intersect the x* = 0 line because of the symmetrical 
nature of the drop. 

Similarly, the radius of curvature Rj can be expressed in terms of x* and cp by 
considering the distance traveled along an elemental length change dS*. Mathematically, 

dS* = Rj d cp (A14) 

However, the arc dS* is perpendicular to the arc traversed by R|; thus, from geomet- 
rical considerations 

dx* = cos cp dS* (A15) 

Combining equations (A14) and (A15) and eliminating dS* yields 


and 


1 1 * 
— + — = z* 

R * n ^ 

1 Xva 


r i 

*t-T 


R 0 

Rp = — 

z X 


z* = - 
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1 deep d sin cp 

Rj dx* dx* 


(A16) 


Substituting the values of Rj and R^ from equations (A13) and (A16) into the 
governing equation (eq. (A10)) yields 


d sin (p f sin cp _ z * 
dx* x* 


(A17) 


However, sin cp can be expressed in terms of the variables z* and x* by noting that 


dz* 

dx* 


= tan cp = 


sin cp 


(A18) 


•V 1 - 


sin 2 <p 


and 


sin (p = 


tan cp 


fc^l + 


(A19) 


- tan (p 


The sign is determined by the quadrant in which the integration occurs. 

Eliminating sin cp from governing equation (A17) by use of equation (A19) and re- 
placing tan cp by dz*/dx* from equation (A18) yield 


d 2 z* JL_ 
dx* 2 x * 


1 + (***\ 

dx*/ 


dz* 

dx* 


= ±z* 


1 + 


dz* 

dx* 


n 3 / 2 


(A20) 


This is the basic governing differential equation for the shape of a drop. The equation is 
a second-order nonlinear differential equation which, up to the present time, has not been 
solved explicitly. 

Bashforth and Adams (ref. 4) solved the same type of equation by a direct numerical 
assault. However, Timoshenko and Woinowsky-Krieger (ref. 8) suggest the introduction 
of the new variable 


j = sin cp 


(A21) 


as a simpler approach to solving the problem. 

Substituting the new variable j (eq. (A21)) into equations (A17) and (A18) yields 
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(A22) 


_dj_ + ± = z * 

dx* x* 


dz* 

dx* 



(A23) 


Thus, the second-order governing equation (A20) has been transformed into two 
first-order differential equations which are much simpler to handle mathematically. 
These equations are equations (6) and (7) in the main text. 
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APPENDIX B 


NUMERICAL PROCEDURES 

The numerical procedures used in solving the governing equations (6) and (7) are de- 
scribed briefly. These equations are solved by a double-precision fourth-order Runge- 
Kutta numerical integration. The numerical procedure for solving two simultaneous 
first-order differential equations by the Runge-Kutta method is given in most texts in nu- 
merical analysis and need not be discussed herein. However, there are a few salient 
points that will be considered. 


Increment Size 

The increment size for the various drop volumes was chosen to be sufficiently small 
such that in most cases the results are accurate to three significant figures. This accu- 
racy was accomplished by iteration using ever smaller values of Ax. The values of Ax 
chosen are listed in table I along with the numerical results. 

To increase the accuracy of the integration in the region where cp equals 90°, Ax 
was reduced by a factor of 100 (0.01 Ax) for 85° < cp < 95°. 

Initial Conditions on z£ 

At the top of the drop, principle radii of curvature become equal because of sym- 
metry. Consequently, at the top of the drop the governing equation (eq. (3)) takes the 
form 

R* = — (Bl) 

° z o* 

where o denotes the apex position of the drop. At the start of the numerical integration, 
the position of the apex of the drop is chosen arbitrarily at some value of z* (fig. 3). By 
choosing z* as an initial condition of the solution, the radius R* and the volume of the 
drop will be uniquely determined. Thus, the drop volume becomes a dependent variable. 

K z* is initially chosen to be small, the radius of curvature at the apex of the drop 
becomes large (see eq. (Bl)). The drop takes on a pancake or disk-like shape. If z* 
is initially chosen to be large, R* becomes small, and the drop takes on a spherical 
shape. By choosing numerous numerical values of z* and solving each problem, the 
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shape of the drop can be found for various values of V (see eq. (A9)). 

In addition, it should be noted that the Runge-Kutta technique puts an upper limit on 
the choice of Ax, in that at the initial starting point, 

Ax* < A- (B2) 

in order to prevent the occurrence of imaginary numbers resulting from the denominator 
in equation (7). 

Thus, for large initial choices of z Q consideration must be given to the choice of 

AX*. 


Tan (p 

The differential dzVcte* in equation (7) equals tan cp. Consequently, for cp < 90° 
the tan cp is positive while for cp > 90°, tan cp is negative. Thus, a ± sign is dis- 
played in equation (7) to indicate that a sign change is required during the numerical inte- 
gration. 


Initial Conditions 

From the choice of coordinates shown in figure 2, the initial conditions for this prob- 
lem are 


X 

* 

II 

o 

z* = 

z o 

(B3) 

= 0 (or j 

= 0) 

^ = 0 

(B4) 

dx* 



i = 

z* 


x* = 0 

0 

(B5) 


x 2 

Initial condition (B5) comes directly from combining equation (A13) with equation (Bl) at 
the apex of the drop. 
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Drop Volume 


Since the solution of the governing equations requires numerical integration, the nu- 
merical results will be a set of ordered numbers z*, x* which describe the shape of the 
drop. The dimensionless volume V* is computed by using the trapezoidal rule as 

v* = l 

6 i=l 

where n represents the total number of numerical steps in the integration. 


+ x i_i + x i 
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APPENDIX C 


GEOMETRICAL RELATION FOR SPHERICAL DROPS 

The volume of a truncated sphere (often called a spherical segment) is given by 

V = — h 2 (3R - h) (Cl) 

3 

(see fig. 10) where h is the vertical distance from the plate to the apex of the spherical 
drop and R is the radius of the sphere. From geometry a unique relation exists among 
R, h, and the contact angle 0; thus, 


R = 


h 

(1 - cos 0) 


Substituting equation (C2) into equation (Cl) and rearranging gives 


or 


h _ h* 

yl/3 v *l/3 



h* 

v *l/3 


3(1 - cos 0) 


ni/3 


7t( 2 + cos 6 ) 


This relation is valid for all contact angles. 

Combining equations (C2) and (C3) gives for the radius of a sphere 


R* 


y*l/3 


1/3 

7 T ' 

(i - cos ef - 1 (l - cos e) 3 

1/3 


3 



(C2) 


(C3) 


(C4) 


(C5) 
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TABLE I. - NUMERICAL VALUES OF DROP SHAPE 
sin 6 y* r* h* h*/V**/^ 

z Q = 0. 4E 02; Ax = 0. IE-05 

0.17460E-01 0. 9 L2b - 1 1 0.873E-03 0.762E-05 

0.8716 Oh -01 0. 568b- 08 0.436E-02 0.190E-03 

0.25883b 00 0.451E-06 0.129E-01 0.170E-02 

0.42262b 00 0.334c-05 0.211E-01 0.468E-02 

0.70712E 00 0.304E-04 0.354E-01 0.14t>E-01 

0.10000E 01 0. 2 6l E - 03 0.500E-01 0.499E-01 

0.70709E 00 0.492b- 03 0.354E-01 0.851E-01 

-0.2 60 43E- 04 U.522E-03 0.205E-02 0.996E-01 

h= 0 . 66 5b - 0 1 h/V* 1 / 3 = 0.826 

z = 0. IE 02; Ax = 0. IE-04 

0. I 7900b- 01 0 . 5 8 9b - 09 0.350E-02 0.306E-04 0.365E-01 

0.8 72U3E- 0 1 0 . 364b - 06 0.1 74E-01 0.762E-03 0.107E 00 

0.2 5884b 00 0.283E-04 0.518E-01 0.681E-02 0.222E 00 

0.42264b 00 0.213E-03 U.844E-01 0.L87E-01 0.313E 00 

0.7G714E 00 0. L 93b - 02 0.14Lt 00 0.583E-01 0.469E 00 

0.10000E 01 0.163E-01 0. 1 99h 00 0.197E 00 0.777E 00 

0.70711b 00 0.303E-01 0.143b 00 0.331E 00 0.106E 01 

— 0. I 0494E- 06 0. 322t-01 0.314E-01 0.382E 00 0.120E 01 

h= 0.260E 00 h/y+i/3 

= 0.816 

z q = 0. 4E 01; Ax = 0. IE-04 

0.1 7460E- 0 1 0 . 9 1 2b - 08 O.8 73E-02 0.762E-04 

0.8 716 lb-0 1 0. 5t>8E - 05 0.436E-01 0.190L-02 

0.25882E 00 0 . 44 8b - 03 0. 129b 00 0.170E-01 

0. 42264E 00 0.328E-02 0.2 1 0b 00 0.464E-01 

0.70711b 00 0.267E-01 0.348b 00 0.143E 00 

0.10000E 01 0.222E 00 0.481b 00 0.458E 00 

0.70709b 00 0.396b 00 U.371E 00 0.734E 00 

-0.2 406 9E- 04 0.426b 00 0.1o7E 00 0.827b 00 

h= 0.53 2E 00 h/V* 1 / 3 = 0.775 


z Q = 0. IE 01; Ax = 0. IE-03 


0. L7453E- 

Cl 

0 . 5 8 3 E 

-06 

0. 349E 

-01 

0.305E 

-03 

0.365E 

-0 1 

0. o71 78E- 

0 1 

0. 359b 

-03 

0. 1 74E 

00 

0.757E 

-02 

0 . 1 0 6E 

00 

0 .2 5886E 

00 

0. 2t>0b 

-01 

0. 5 02 b 

00 

0 .650E 

-01 

0.219b 

00 

0. 4?2o5£ 

0 0 

0. 16 5E 

00 

0. 7d2b 

00 

0 . 1 6 6 F 

00 

0 .304E 

00 

0. 70711b 

00 

0. 1 04E 

01 

0. 1 lbb 

01 

0 • 4 39 E 

00 

0.433b 

00 

0. 1 0000b 

0 1 

0. 484E 

01 

U. 167b 

01 

0 . 1 08 E 

01 

0 .636E 

00 

0. 7071Gb 

0 0 

0 . 755b 

01 

0.130b 

01 

0.15 IE 

01 

0.767b 

00 

-0.2 1354b- 

03 

0. 81 6b 

01 

0. 992b 

00 

0.164E 

01 

0 .8 15E 

00 

h = 

0 . 12 IE 

01 


h/V* 1 / 3 = 

0.599 





0.365E-01 
0 . L0 7E 00 
0.222E 00 
0.313E 00 
0 . 466b 00 
0.7 5 7E 00 
0.100b 01 
0.110E 01 


0 .365E-0 1 
0.107E 00 
0.222E 00 
0.313E 00 
0.469E 00 
0.781E 00 
0.108E 01 
0.124E 01 
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TABLE I. - Continued. NUMERICAL VALUES OF DROP SHAPE 
sin 9 V* r* h* hVV* 1 / 3 


z Q = 0. 4E 00; Ax = 0. IE -03 


0. 1 74 3 7E- 

0 1 

0. 909b- 

-03 

O. 8 72b- 

01 0.761E- 

-03 

0 .365E 

-01 

0.87171E- 

0 1 

0. 52 7t ■ 

-02 

0. 426E 

00 0.184E- 

-01 

0 . 1 06E 

00 

0. 258856 

00 

0.2 7LE 

00 

0.111b 

01 0.136b 

03 

0 .210b 

00 

0.42263b 

0 0 

0. 123b 

Ol 

0.157b 

01 0 . 2 9 9 E 

00 

0 .2796 

00 

0. 7G711E 

uo 

0. 5 03b 

01 

0.208b 

Ol 0.549E 

00 

0 .3786 

00 

0. i OOOOt 

0 L 

0 . 1 6 7b 

02 

0. 24 lb 

01 0 . 1 3 3 E 

01 

0.527E 

00 

0. 707 10b 

0 0 

0 . 2 4 5b 

02 

0.2 2*3 b 

01 0.179E 

01 

0.618b 

00 

-0.227436- 

0 3 

0 • 2 o 3b 

02 

0. 190E 

01 0.194b 

01 

0.650E 

00 

h = 

U. 146b 

0 l 


h/V* 1 / 3 = 0.438 







Z Q = 0. IE 

00; Ax 

= 0. IE-03 




0. L 7456b- 

01 

0.553b 

-03 

0. 344E 

00 0.298E- 

-02 

0 .36 3E 

-0 l 

0.871366- 

ol 

0 . 1 69t 

00 

0. 138E 

01 0.540b- 

-01 

0.977b 

-01 

0. 2 5a a 3E 

00 

0.2 6 9E 

01 

0. 253b 

01 0.239E 

00 

3.172E 

00 

0.42262b 

0 0 

0. 743t 

01 

0.3 06b 

01 0.429E 

00 

0.220E 

00 

0. 707 12b 

0 0 

0. 205b 

02 

0. 3 6 Ob 

01 0.7956 

00 

0.291E 

00 

0.1 00 OOF 

0 l 

0. 32ob 

02 

0.393b 

01 0.149E 

01 

0 . 399b 

00 

0.707026 

0 0 

0. 736b 

02 

0. 375b 

01 0.1946 

01 

0 . 463b 

00 

-0. 19288b- 

03 

0.798b 

02 

0. 341b 

01 0.209b 

01 

0 .485E 

00 

h = 

0 . 1 6 3E 

01 


h/V* 1 / 3 = 0.382 







Z Q = 0. 4E 

1 

o 

h-i 

& 

= 0. IE-03 




0.1 7432b- 

01 

0 . 6 9oE 

- 02 

0.8 06 E 

00 0.67&E- 

-02 

0.354E 

-01 

0.87138F- 

Cl 

0.73 Ob 

00 

u. 2 34b 

01 0.768E 

-01 

0.853b 

-01 

0.2 5884b 

00 

0 . 642b 

01 

0.334b 

01 0.269E 

00 

0 . 145b 

00 

0 . 4 ?2o 4b 

0 0 

0. 1 52b 

02 

0.407b 

01 0.458E 

00 

0 . 185E 

00 

0. 70717b 

00 

0. 3 7 Ob 

02 

0.46 lb 

01 0.821E 

00 

0.246E 

00 

0. 100 00b 

01 

0. 873b 

02 

0. 493b 

01 0 . 1 51E 

01 

0.340E 

00 

0.70704b 

Oo 

0. 12lb 

03 

0.476b 

01 0.196E 

01 

0 . 396E 

00 

-0. 1 2740b- 

0 3 

0.131b 

03 

0.441b 

01 0.21 IE 

01 

0 .415E 

00 

h = 

0 . 1 7 lb 

01 


h/V* 1 / 3 = 0.337 







Z Q = 0. IE 

-01; Ax 

= 0. IE -02 




0.1 7455b- 

01 

0. 1 09b 

00 

0. 2 1 Ob 

01 0.145E- 

-01 

0 .3036 

-o l 

0.8 7176b- 

0 l 

0.267E 

01 

0. 3b7b 

01 0.916E- 

-01 

0 . 660E 

-01 

0.2 5886b 

00 

0. 1 52b 

02 

0.5 06E 

01 0.282E 

00 

0.114E 

00 

0.42286b 

00 

0. 3 2 0b 

02 

0.538b 

01 0.468E 

00 

0 . 147E 

00 

0. 70713b 

00 

0. 707b 

02 

0.61 lb 

01 0.824E 

00 

0 . 199E 

00 

0.999996 

oo 

0. 136E 

03 

0. 643b. 

01 0.150E 

01 

0.279E 

00 

0 . 7 03 3 66 

00 

0.2 i3E 

0 3 

0.625L 

01 0.195E 

01 

0.326E 

00 

-0.6 796 3E- 

03 

0. 230E 

03 

0. 590E 

01 0.210E 

01 

0.342E 

00 

h = 

0 . 177E 

01 


h/V* 1 / 3 = 0.2 89 
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TABLE I. - Continued 


NUMERICAL VALUES OF DROP SHAPE 


sin d V* r* h* 


z Q = 0. 4E-02; Ax = 0. IE-02 


0. 1 7462b- 

-01 

0.3 08b 

00 


0.31 IE 

01 

0.1 74E 

-01 

0.87228b- 

-01 

0. 466b 

01 


0. 48 BE 

01 

0.940 E 

-01 

0.25898b 

00 

0.231b 

02 


0.605E 

01 

0.282E 

00 

0.42270E 

00 

0.463b 

02 


0.657E 

01 

0.466E 

00 

0. 7 0 7 l 6£ 

00 

0. 986E 

02 


0. 709b 

01 

0.8 19E 

00 

0.99999E 

00 

0.212E 

03 


0.74 IE 

01 

0. 150E 

01 

0.70676b 

00 

0.28 7E 

03 


0. 723b 

0L 

0. 194E 

01 

— 0.4877 5b- 

-03 

0. 31 lb 

03 


0.688b 

01 

0 .209E 

01 

h = 

o 

• 

p~ 

c 

m 

01 



h/V* 1 / 3 = 

= 0.256 




N 

O 

II 

0. IE 

-02; Ax 

= 0 . 

IE-02 


0 . L 7460E- 

•01 

0.826b 

00 


0 . 464b 

01 

0. 187E- 

-01 

0.8 7182E- 

-01 

0. 86 5b 

01 


0.638b 

01 

0.941E- 

-01 

0.25890b 

00 

0.377b 

02 


0. 753b 

01 

0.279E 

00 

0.42290E 

00 

0. 72 7b 

02 


0. 805b 

01 

0.462E 

00 

0. 707 53E 

00 

0. 149c 

03 


0. 856b 

01 

0.8 12E 

00 

0. 10000b 

0 l 

0.312b 

03 


0.888b 

01 

0. 148E 

01 

0.7058 IE 

00 

0.41 9b 

03 


0.870b 

01 

0. 192E 

01 

-0. 1 40 43b- 

02 

0.454b 

03 


0. 835b 

01 

0.207E 

01 

h = 

0 . 1 a 4E 

01 



h/V * 1 / 3 = 

= 0.239 




Z 

'o " 

0. IE 

o 

CO 

& 

= 0 . 

IE-02 


0.1 7463b- 

0 1 

0.222L 

01 


0. 712b 

01 

0. 18 7 E- 

-01 

O.d 7192E- 

01 

0. 178b 

02 


0.882b 

01 

0.926b- 

-01 

0.25901E 

00 

0.697b 

02 


0. 996b 

01 

0.275E 

00 

0.42281b 

00 

0. 1 2 9E 

03 


0 . 1 U5b 

02 

0.455E 

00 

0. 70778b 

00 

0.255b 

03 


0.110b 

02 

0.802b 

00 

0.99999E 

00 

0.51 7 E 

03 


0.113b 

02 

0.147b 

01 

0.7 053 0b 

00 

0. 691b 

03 


0. 1 L lc 

02 

0.191E 

01 

-0.71106b- 

0 3 

0. 748b 

03 


0.108b 

02 

0.206E 

01 

h = 

0. 18 7E 

01 



h/v * 1 / 3 - 

; 0.206 






IS) 

o 

II 

© 

IE-05; Ax 

= 0 . 

IE -02 


0. L 7457E- 

01 

0.688b 

01 

0. 120b 

02 

0.182b 

-01 

0.87192E- 

01 

0.454b 

02 

0. 1 36b 

02 

0.907E 

-01 

0.25898b 

00 

0. 1 6 1 b 

03 

0. 1 48b 

02 

0.271 E 

00 

0.42270b 

00 

C 

• 

00 

-si 

T 

03 

0. 1 52b 

02 

0.448b 

00 

0.70768b 

00 

0. 547b 

03 

0. 158E 

02 

0.791 E 

00 

0. 10000E 

01 

0. 1 U8t 

04 

0. 1 6 IE 

02 

0.145b 

01 

0. 70600b 

00 

0. 143E 

04 

0. 1 5 9b 

02 

0.189b 

01 

-0. 1 4109E- 

02 

0. 1 55b 

04 

0.155b 

02 

0.204E 

01 


h = 0.19 it oi h/V* 1 / 3 = 0.16 5 


h^yv * 1 / 3 


0 .257E-01 
0.563E-01 
0.990E-01 
0.130E 00 
0.177E 00 
0.251E 00 
0.294E 00 
0.308E 00 


0 .200E-0 l 
0. 458E-01 
0.83 3 £-01 
O.lllE 00 
0.153E 00 
0.219E 00 
U.257E 00 
0.270E 00 


0. 144E-01 
0 .354E-01 
0 . 669t-0 1 
0.901E-01 
0.126E 00 
0.183E 00 
0.2L6E 00 
0.227E 00 


0 .959E-02 
0.254E-01 
0.497b-01 
0.679E-0L 
0.967E-01 
0.142E 00 
0.L68E 00 
0.176E 00 
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TABLE I. - Concluded. NUMERICAL VALUES OF DROP SHAPE 


sin 0 


h* 


h*/V 


*1/3 


z Q = 0. IE-07; Ax = 0. IE-02 


0. 1 74556- 

■01 

0. 140E 

02 


0. 1 o 76 

02 

0. 130 E 

-01 

0 .747E-02 

0.8 71926- 

01 

0.851E 

02 


0. 184E 

02 

0.893 E 

-01 

0 .2046-01 

0. 2 59 0 IE 

00 

0 . 2 8 8t 

03 


0. 1 95t 

02 

0 . 2 o 8 E 

00 

0 .406E-01 

0.422996 

00 

0.5 G4E 

03 


0. 200fc 

02 

0 . 444 E 

00 

0.558E-01 

0. 70769E 

00 

0. 944E 

03 


0.2056 

02 

0.785E 

00 

0 .8006-01 

0.99999E 

00 

0. 183E 

04 


0.2 086 

02 

0.144E 

01 

0.1186 00 

0. 7 060 'IE 

00 

0.242E 

04 


0.2 006 

02 

0.183E 

01 

0.140E 00 

-0.8 689 2£- 

03 

0 . 2 1> 2 1 

04 


0. 2 036 

02 

0 .2036 

01 

0.1476 00 

h = 

0 . 1 9 36 

01 



h/v* 1 / 3 = 

0.140 






N 

O 

II 

0. IE 

1 

o 

CD 
>* • 

& 

= 0. 

IE -02 



0. 1 7466E- 

0 l 

0 . 2 3 d 6 

02 


0.2 146 

02 

0.179E 

-01 

0 .625E-02 

0.8 71 94F- 

0 1 

0. 137E 

0:5 


0.2316 

02 

0 .B92E- 

-01 

0 . 173E-0 l 

0 . 2 58 88h 

uo 

0 . 4 5 1 E 

03 


U.2426 

02 

0 .267E 

00 

0 .348E-01 

0.422886 

00 

0. 7 8 1 1 

03 


0.2476 

02 

0.4426 

00 

0.480E-01 

0. 707556 

00 

0. 1 456 

04 


0.2526 

02 

0. 78 l c 

00 

0.691E-01 

0. 99949E 

00 

0.2 77E 

04 


0.2556 

02 

0. 144-E 

01 

0.102E 00 

0. 705 80E 

0 0 

0. 36oE 

04 


0.253E 

02 

0.187E 

01 

0.121E 00 

-0. 198 54E- 

02 

U. 3 96E 

04 


0.2506 

02 

0.202E 

01 

0.128E 00 

h = 

0 . 1 9 4E 

01 



h/V* 1 / 3 = 

0.123 






N 

O 

II 

0. IE 

-12; Ax 

= 0 . 

IE-02 



0. 1 74626- 

0 1 

0.422E 

02 


0.285E 

02 

0. 178E 

-01 

0. 51 IE-02 

0. 87209E- 

01 

0. 2 37b 

03 


0. 3 ulE 

02 

0 .888E- 

-01 

0 . 144E-0 1 

0.25895E 

0 0 

0. /t>3h 

03 


0. 3126 

02 

0.265E 

00 

0.291E-01 

0.422636 

0 0 

0. 1 3 lb 

04 


0.31 76 

02 

0.440E 

00 

0 .402E-01 

0 . 707836 

0 0 

0.239E 

04 


0.3226 

02 

0.773E 

00 

0 .582E-01 

0.99999E 

00 

0. 455E 

04 


U. 32 56 

02 

0.143E 

01 

0 .863E-0 l 

0.70552E 

00 

0 . 5 99E 

04 


0.3236 

02 

0. 186E 

01 

0.1036 00 

-0. 1 2490c- 

0 2 

U. b48E 

04 


0. 3206 

02 

0.202E 

01 

0.108E 00 


h= 0.19 5E 01 


h/V * 1 / 3 = 


0.105 
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C-66-3242 


(c) 1- g's. 

2 

Figure 1. - Variation of mercury droplet 
shape with gravity level. (Ref. 3) 





(a) Wetting drop, O°<0<9O°. 


(b) Intermediate case, 0 ■ 90°. 


(c) Nonwetting drop, 
90° < 0 < 180°. 


(d) No surface attraction. 


Figure 2. - Drops resting on flat surface, shown in order of decreasing molecular attraction between liquid and solid. 


Reference plane 







Dimensionless volume, 






Figure 9. - Relation between contact and maximum drop 
radius. 
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